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Results of investigation of radio galaxies of the survey 
"Cold": photometry, colour redshifts and the age of the 
stellar population 
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Abstract. BVRI data for the majority of the objects of the RC catalogue with steep spectra 
and m,R < 24 m are presented. These data have been used to estimate colour redshifts and the age 
of stellar systems of host galaxies. By way of example of distant radio galaxies it is shown that 
this approach gives an accuracy of redshift estimates close to that for field galaxies (~ 20 %). The 
age estimates are less confident, however, the lower age limit for not too distant (z < 1.5) objects 
is determined quite reliably. Several galaxies have been detected that have an age above that of 
the Universe at the given z in a simple CDM model of the Universe. A possibility of using such 
objects to specify the part played by "dark energy" is discussed. This paradox disappears in a 
model with the A-term equal to 0.6-0.8. 
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1. Introduction 

A population of powerful radio galaxies is observ- 
able with the available facilities practically to any 
distances. This allows one to follow the evolution of 
this population in the radio range from the moment 
it originated up to the present time. 

However, a good deal of effort has to be under- 
taken to find new objects at z > 3 with application of 
rigorous methods of selection by radio, optical and in- 
frared properties, which leads to inhomogeneous sam- 
ples with regard to different parameters. 

It is customary to assume that a population of 
very powerful radio galaxies of type FRII owes its 
origin to giant elliptical galaxies having supermassive 
(10 9 Mq) black holes at the centre. That is why the 
evolution of this population is related to the prob- 
lem of formation of the largest stellar systems and 
to the problem of evolution of massive black holes. 
Besides, these objects are often associated with clus- 
ters and groups of galaxies and can be indicators of 
distant clusters formed in the nodes of the large-scale 
structure. At last, there is a number of suggestions for 
using them in estimation of parameters of their en- 
vironment and even of the geometry of the Universe 
and its dynamics (Parijskij et al., 1998; Daly, 1994). 



In contrast to quasar, in radio galaxies one can 
study in details the stellar population. As a rule, 
in objects with medium and weak energy release of 
the active nucleus in the optical and near infrared 
ranges, radiation of stars rather than gas dominates 
at z < 1 — 1.5. For this reason, one may attempt to 
apply the method of stellar evolution and evolution of 
synthetic colours of the stellar population to determi- 
nation of "colour" redshifts and the age of the stellar 
population. When lucky, one can define the moment 
of star formation (z). 

In the generally accepted scheme of the Unified 
Theory of active objects, radio galaxies of type FRII 
differ from quasars, BLLac objects and other popu- 
lations with powerful nucleus activity only in that at 
what angle the axis of rotation of the gas-dust tore, 
that screens the accretion disk around the black hole, 
is located with respect to the observer. And the data 
on evolution of these objects can be applied to other 
classes of objects. 

At last, powerful galaxies with active nuclei (radio 
galaxies and quasars, predominantly "radio quiet") 
supply much UV radiation to the Universe, and in 
combination with the integral optical depth by Thom- 
son scattering of relic photons of the 3K background 
permit the moment of secondary ionization of the 
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Universe to be refined. 

Deep optical investigations of host galaxies in this 
class of objects are impeded by the fact that their spa- 
tial density is by 5-6 orders of magnitude less than 
that of background galaxies, and this is why they are 
virtually lacking in ultimately deep frames of small 
size. In the HDF (Hubble Deep Fields) there is found 
but one distant (z — 4.42) radio galaxy of medium 
power. For this reason, sampling of fields for such 
investigations has to be started with preliminary se- 
lection on radio astronomy data with allowance made 
for all possible indirect criteria. The principles of se- 
lection of candidates for distant objects in the "Big 
Trio" project (Parijskij et al., 1994, 1096a, b) are sim- 
ilar to universally adopted (McCarthy, 1993). The list 
of objects presented below has been derived from the 
catalogue of the "Cold" survey (Parijskij et al., 1991, 
1992) as a result of several selection steps: selection 
of comparatively faint sources having steep spectra 
(a < —1.05 ~ v a ) and being of the FRII morpho- 
logical type (Fanaroff, Riley, 1974) in the cm range 
(10-100 mJy). Further selection was done according 
to the procedure described below. 

The labour input for acquisition of quality spec- 
troscopic data on distant and faint galaxies and radio 
galaxies forces to search for indirect methods for de- 
termination of the redshift and other characteristic 
features of these objects. With regard to the power- 
ful radio galaxies, even photometric estimates proved 
to be helpful and they have so far been widely used 
(McCarthy, 1993; Benn et al., 1989; Parijskij et al., 
2000a, b). 

Over the last few years, colour characteristics of 
faint galaxies have come to be used in addition to 
the methods of photometric evaluation of the redshift, 
and this approach forms a basis for a number of major 
projects (e.g. see Szalay, 1996). As we have already 
noted, the colours of stellar systems make it possi- 
ble to estimate their age as well. Generally speaking, 
the legitimacy of using colours for such estimations of 
evolution characteristics of the population of powerful 
radio galaxies requires a separate study because the 
effect the nuclear activity in these objects has on the 
colour characteristics is not sufficiently understood. 
The sites of "secondary" star formation may also tell 
upon the colour characteristics. 

The present paper is a continuation of our pro- 
gramme of investigation of the radio galaxies detected 
in the RATAN-600 survey "Cold" with involvement 
of multicolour photometry data for assessment of 
colour redshifts and ages of stellar systems of host 
galaxies. Current use of these data is practically in- 
evitable. Direct spectroscopy of faint objects required 
until recently a great deal of the 6 m telescope ob- 
servational time when working with objects fainter 
than 20™. So, in 1995-1996 we managed to measure 
the spectroscopic redshift but in four bright objects 



(three quasars and one galaxy, all being brighter than 
tur = 20. 5 m ) and only in one faint galaxy, jur = 23™ 
(Dodonov et al., 1999). The redshift of the latter, 
z = 2.73, was estimated from the only emission line. 
It is independently confirmed by our colour data: the 
negative B-V colour index agrees with identification 
of a strong line as Ly a , which is typical of radio galax- 
ies. 

During the last time, in connection with placing 
in service of more efficient spectral tools in SAO 
RAS, spectroscopy of faint objects is expected 
to substantionally improve. The spectra of an- 
other two objects, radio galaxy RC J 0908+0451 
and quasar RCJ 1154+0431, were obtained at the 
6 m telescope in 2001 March with the aid of new 
equipment SCORPIO (Afanasiev and Moiseev, 
http://www.sao.ru/~moisav/scorpio/scorpio.htmfy . 
The spectral redshift for the radio galaxy is practi- 
cally the same as the colour redshift. The data on 
the colour redshift will help improve the efficiency of 
spectroscopic observations. 

This paper has the following structure. In the first 
section we present the data of multicolour BVRI pho- 
tometry and also the latest results of identifications of 
the "Cold" survey radio sources. Further we describe 
the procedures of obtaining colour redshifts and ages 
of stellar systems from BVRI magnitudes with the 
use of present-day models of spectral energy distribu- 
tion (SED). In the sections that follow we discuss the 
results obtained, evaluate model differences in the es- 
timates of redshifts and propose new steps to improve 
the accuracy of estimations. 

2. Multicolour photometry of RC ob- 
jects 

By the present time observations of about 60 RC ob- 
jects (radio galaxies and quasars) of our sample have 
been carried out with the 6 m telescope of SAO RAS 
in four filters (B, V, R, I). In this paper we present 
the photometry results for 50 radio galaxies. The ob- 
servations were made in 1994-1998. A CCD ISD015A 
of 580 x 520 pixels with a pixel size of 0.205" x 0.154" 
was employed in the 1994-1995 observations, and a 
CCD ISD017A (1160x 1040, 0.137"x0. 137") was used 
in observations of 1996-1998. When reading the lat- 
ter, the initial elements of the matrix were summed 
up to produce the output elements twice as large in 
both coordinates. The exposure time was defined in 
the course of observations, proceeding from bright- 
ness and colour of objects, from 1200-1800s (total) 
in the B band to 400-800 s in the R band for a typical 
object of our sample having R = 22 m — 23 m , so that 
the resulting signal-to-noise ratio would be no worse 
than 4-5 in all the filters. At a seeing worse than 2", 
we had to increase the exposure time by a factor of 
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Tabic 1: A list of RC objects classified as QSO and 
excluded from this study 



Name RA(2000.0) Dcc(2000.0) 



RC 


J0038+0449 


00 


38 


34.65 


+04 


50 50.5 


RC 


J0042+0504 


00 


42 


27.17 


+05 


05 24.7 


RC 


J0126+0502 


01 


2G 


16.11 


+05 


02 09.9 


RC 


J0143+0505 


01 


43 


33.89 


+05 


07 57.4 


RC 


J0226+0512 


02 


2G 


19.80 


+04 


46 32.5 


RC 


J0459+0456 


04 


•59 


04.28 


+04 


55 54.4 


RC 


J0506+0558 


05 


0G 


25.00 


+05 


08 19.3 


RC 


J1154+0431 


11 


54 


53.50 


+04 


24 12.5 


RC 


J1740+0502 


17 


40 


33.96 


+05 


02 42.3 


RC 


J2013+0508 


20 


13 


23.48 


+05 


10 30.5 


RC 


J2036+0451 


20 


3G 


56.93 


+04 


49 52.7 



1.5-2. 

Identifications of RC objects are presented in 
Fig. 5 (see after the references). The radio isophotes 
were plotted from the data of VLA observations and 
superimposed on the images of galaxies obtained with 
the 6 m telescope in the R filter. The VLA radio ob- 
servations were carried out in the 1990s with differ- 
ent antenna configurations (Parijskij et al., 1996a). 
Besides, the VLA data for the RC objects from the 
MIT survey (Fletcher et al., 1996) and the FIRST 
survey (White et al., 1997) data were used. Images 
are also presented for the objects that we classified as 
QSO and which are excluded from the present study. 
The list of the objects is given in Table 1. 

The processing of optical images was accom- 
plished using the standard procedure in the MIDAS 
system. Subtraction of the averaged dark frame and 
element-by-element correction with the use of twilight 
sky frames were performed. The residual background 
inhomogencity in the I filter, which is due to inter- 
ference of night sky emission lines, was removed by 
means of the procedure of subtraction of the median 
sum of all working frames of a given night or sev- 
eral nights of one observing run. The calibration of 
photometric measurements for their conversion to the 
standard Johnson-Cousins system was done using the 
stars from Landolt's (1992) list, which were observed 
several times during a night. To do photometry of 
the selected object, a circular aperture of one and the 
same size was used when measuring in different filters. 
The aperture size was chosen to be between 3" and 
12", depending on the object luminosity. The typical 
size was 4" — 5". The background was measured with 
a circular aperture of sufficiently large radius in or- 
der not to cover the outer regions of the object being 
measured. When necessary, the neighbouring objects 
were removed by the procedure of interpolation of 
the surrounding background. The photometry accu- 
racy was generally not worse than 0.1 m for galaxies 
brighter than 21 m . It dropped down to 0.2 m - 0.25 m 



Table 2: Characteristics of photometric bands 



Filter name 


A eJ7 


A/E(B - V) 


C 


Landolt B 


4400 


4.315 


3.620 


Landolt V 


5500 


3.315 


3.564 


Landolt R 


6500 


2.673 


3.487 


Landolt I 


8000 


1.940 


3.388 



for 23 m - 24 m and was as low as 0.3 m - 0.5 m for 
galaxies fainter than 25 m — 25. 5 m . 

To provide for the Galactic extinction, we used 
the charts from the paper "Maps of Dust IR Emis- 
sion for Use in Estimation of Reddening and CMBR 
Foregrounds" (Schlcgel et al., 1998), written in the 
form of FITS files. The conversion of stellar magni- 
tudes to flux densities was implemented via the for- 
mula S(Jy) = l0 c -°- 4m (von Hocrner, 1974). The 
values of the constant C for different filters are listed 
in Table 2, where the following characteristics are also 
given: name of filter, effective wavelength, coefficient 
A/E(B — V) of change-over from distribution of dust 
radiation to extinction in a given band under the as- 
sumption of extinction curve Ry = 3.1. 

The stellar magnitudes of 50 galaxies of our sam- 
ple, which were corrected for extinction in the four 
filters, are presented in Table 3. 

3. Models of energy distribution in the 
spectra of host galaxies 

In the late 1980s and early 1990s, attempts were made 
to use colour characteristics of radio galaxies for esti- 
mation of rcdshifts and ages of systems of host galax- 
ies. There appeared numerous evolution models with 
which observational data were compared and results 
largely differing from one another were obtained (Ari- 
moto & Yoshii, 1987; Chambers & Chariot, 1990; 
Lilly, 1987, 1990; Parijskij et al., 1996a). Over the past 
few years the models like PEGASE: Project de'Etudc 
des Galaxies par Synthcse Evolutive (Fioc & Rocca- 
Volmerange, 1997), and GISSEL'98 (Bruzual, Char- 
lot, 1993; Bolzonella et al., 2000), in which an attempt 
was made to eliminate the shortcomings of the pre- 
ceding versions, have been most commonly used. 

In the "Big Trio" experiment (Parijskij et al., 
1996a), we also tried to apply these methods to dis- 
tant steep spectrum objects of the RC catalogue, for 
which we measured the magnitudes in the four filters 
(BVRI). This paper is distinguished from the paper 
of Parijskij et al. (1996a) by using more reliable iden- 
tifications of radio sources, new photometric data and 
new models of energy distribution. Besides, we used 
in our work the smoothing procedure that makes it 
possible to model and predict the flux in the given 
filter at the specified SED with allowance for the fil- 
ter throughput curve and also for the redshift effect. 
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Table 3: Stellar magnitudes of the sample radio galaxies corrected for extinction 



Source 


B 


V 


R 


I 


Source 


B 


V 


R 


I 




m 


m 


m 


m 




m 


m 


m 


m 


0015+0503a 


23.89 


22.97 


22.20 


21.36 


1155+0444 


21 


36 


19.83 


18.90 


18.20 


0015+0501 


24.82 


23.91 


23.37 


22.22 


1213+0500 


23 


55 


22.90 


22.04 


21.32 


0034+0513 


25.28 


24.79 


23.25 


21.79 


1235+0435b 


24 


15 


22.81 


21.59 


20.35 


0039+0454 


24.81 


24.00 


22.69 


21.22 


1322+0449 


23 


68 


22.52 


20.77 


19.17 


0105+0501 


24.00 


22.48 


22.78 


22.43 


1333+0452 


24 


87 


24.44 


23.56 


22.46 


0135+0450 


20.49 


19.16 


18.42 


17.82 


1339+0445 


25 


05 


23.72 


22.70 


21.55 


0152+0453 


23.31 


23.02 


22.47 


21.70 


1357+0453 


22 


98 


21.85 


21.10 


20.11 


0159+0448 


22.65 


21.72 


21.23 


20.65 


1429+0501 


25 


57 


23.24 


21.64 


20.50 


0209+0501a 


20.37 


19.19 


18.43 


17.78 


1436+0501 


23 


90 


23.86 


23.39 


22.70 


0209+0501b 


25.72 


24.09 


23.12 


21.63 


1446+0507 


21 


48 


20.03 


19.17 


18.54 


0318+0456 


25.64 


23.98 


22.61 


20.99 


1503+0456 


24 


02 


23.67 


23.14 


22.24 


0444+0501 


23.48 


23.70 


23.33 


23.13 


1510+0438 


24 


98 


23.73 


22.57 


21.25 


0457+0452 


22.01 


20.86 


20.05 


19.37 


1551+0458 


25 


57 


25.34 


24.43 


23.30 


0836+0511 


23.68 


23.53 


23.09 


22.44 


1626+0448 


22 


32 


23.07 


22.73 


22.63 


A007 i C\ A AC 1 

0837+0446 


no f\o 

23.03 


no on 

23.29 


oo nn 

22.99 


oo 1 1 

22.11 


1638+0450 


22 


86 


oo oo 
22.66 


22.14 


o 1 n a 

21.04 


0845+0444 


24.72 


22.42 


21.09 


19.77 


1646+0501 


24 


01 


22.44 


20.97 


19.76 


0908+0451 


21.63 


20.72 


19.85 


19.07 


1703+0502 


24 


22 


23.39 


23.12 


22.26 


0909+0445 


22.60 


21.53 


20.50 


19.59 


1706+0502 


24 


73 


24.19 


23.25 


21.88 


0934+0505 


25.29 


24.45 


24.67 


23.61 


1722+0442 


22 


30 


21.59 


20.63 


19.44 


1011+0502 


23.71 


23.18 


22.47 


22.60 


2029+0456 


22 


85 


22.24 


21.66 


20.53 


1031+0443 


23.93 


22.79 


22.09 


20.85 


2219+0458 


24 


80 


25.03 


23.72 


22.25 


1043+0443 


23.98 


23.57 


22.51 


21.70 


2224+0513 


23 


16 


22.31 


21.43 


20.32 


1124+0456 


20.30 


18.79 


17.85 


17.07 


2247+0507 


23 


64 


23.18 


22.53 


21.43 


1142+0455 


24.83 


22.53 


21.38 


20.39 


2348+0507 


23 


89 


23.79 


23.56 


23.08 


1152+0449 


23.86 


23.66 


22.39 


21.00 
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Figure 1: Hubble diagram R-value vrs redshift for different radio galaxies taken from literature. The figure is 
reproduced from Pursimo et al. (1999). 
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These changes in the procedure enabled making the 
results of estimation more reliable as compared to our 
previous work. 

We (Verkhodanov et a!., 1999) preliminarily dis- 
cussed applicability of the new methods to a popu- 
lation of distant (z > 1) radio galaxies with known 
redshifts, for which we managed to find in literature 
more or less reliable data of multicolour photometry 
in the optical and near infrared ranges in no fewer 
than three filters. It is shown, in particular, that red- 
shifts can be estimated with an accuracy of 25-30 % 
at 1 < z < 4, given measured stellar magnitudes in 
more than three filters. And given at least one lu- 
minosity estimate in the infrared region, it suffices 
to use measurement in three filters. Estimates were 
made for two evolution models. 

As one of the SED models the evolution model 
PEGASE (Fioc & Rocca-Volmerange, 1997) was em- 
ployed for galaxies of the Hubble sequence both star 
forming and evolving in a passive fashion. One of the 
principal merits of this model is the extension to the 
near infrared range (NIR) of the atlas of synthetic 
spectra of Rocca-Volmerange & Guiderdoni (1988) 
with the revised stellar library including parameters 
of cool stars. The model covers a range from 220 A to 
5 microns. According to the authors the algorithm of 
the model makes it possible to follow rapid evolution 
phases such as red supergiants or AGB in the near IR. 
We used a wide set of SED curves from this model in 
a range of ages from 7 x 10 6 years to 19 x 10 9 years 
for massive elliptical galaxies. 

From the library of synthetic spectra of the model 
GISSEL'98 (Bolzonella et al., 2000) derived with the 
aid of the evolution models of Bruzual & Chariot 
(1993, 1996), we used the computations for elliptical 
galaxies. The library of synthetic spectra was con- 
structed with the following star formation parame- 
ters: simple stellar population (SSP), the duration of 
starburst activity is 1 billion years, the starburst ac- 
tivity decays exponentially. The solar metallicity was 
used in the model. The initial mass function (IMF) 
with an upper limit of 125 solar masses was taken 
from Miller & Scalo (1979). As Bolzonella et al. (2000) 
show the choice of IMF has no effect on the accuracy 
of redshift measurements. The mode tracks are com- 
puted in a wavelength rang e from 200 to 95800 A. We 
used the range specified by a redshift limit of to 6 
in our calculations. 

The sets of evolution models employed are avail- 
able at http://sed.sao.ri\ (Verkhodanov et al., 2000). 



4. Estimation of age and redshift 

Before using the model transmission curves we per- 
formed their smoothing with the filter via applying 



the algorithm 

n 

Zj s i-n/2+jfjk{z) 

a _ J^o 

°ik — n i 

3=0 

where Si is the initial model curve of SED, Sik is 
the model curve of SED smoothed by the fcth filter, 
fk(z) is the curve of transmission of the fcth filter 
"compressed" (1 + z) times when "moving" along the 
wavelength axis of the SED curve, j = 1, n is the num- 
ber of the pixel in the filter transmission curve. From 
the fc curves of SED thus formed (there are four of 
them in our case), we constructed a two-dimensional 
array (A, filter) of smoothed spectra for further com- 
putations. 

The evaluation of ages and redshifts was carried 
out by the method of choosing an "optimum" location 
of photometric magnitudes obtained in photometric 
observations of radio galaxies in different filters on 
the SED curves. We used the already computed and 
table-specified SED curves for different ages. The al- 
gorithm of selection of the optimum location on the 
curve was briefly as follows (for details see Verkho- 
danov, 1996): by way of shifting the points lengthwise 
(along the wavelength axis) and transverse (along the 
intensity axis) the SED curves, we defined the lo- 
cation at which the sum of squares of deviations of 
points from the corresponding smoothed curves is a 
minimum, i.e. the minimum of x 2 is actually calcu- 
lated: 



N filters 

E 

k=l 



F, 



bs.k 



p-SED fc (z) 



where F b St k is the observed stellar magnitude in the 
fcth filter, SEDfe(z) is the model stellar magnitude for 
the given spectral distribution in the fcth filter at a 
given z, p is the free coefficient, ak is the measurement 
error. 

The redshift was found from the shift of location of 
the observed magnitudes at their best location on the 
SED curves from the "rest frame" position (position 
at z = 0). From the total set of curves for different 
ages we chose the ones on which the sum of squares 
of deviations for the given observations turned out to 
be a minimum. 

Fig. 6 (see at the end of the paper) gives the results 
of the most plausible choice of evolution models for 
the photometry data and the corresponding function 
of plausibility. 

5. Photometric redshifts for RC cata- 
logue radio galaxies 

To examine the potentialities of the method for deter- 
mination of redshift and ages of the stellar population 
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of host galaxies from the data of multicolour photom- 
etry, we selected 40 remote radio galaxies with known 
redshifts, for which stellar magnitudes in no fewer 
than three filters are available in literature (Verkho- 
danov et al., 1998b, 1999). 

It should be noted that literature photometric 
data are most inhomogeneous: they were obtained 
by different authors, different tools, with different fil- 
ters; measurements for one and the same object were 
not always made in the same apertures, etc. This is 
why, only 42 out of 300 radio galaxies of the original 
sample remained. The subsamples of objects (3C, 4C, 
B2) with relatively homogeneous data are not large 
enough to be statistically compared. At first, using 
the collected photometric data obtained with the aid 
of the models PEGASE and GISSEL'98, only the ages 
of the stellar population of host galaxies with a fixed 
known rcdshift were derived. After that, a search for 
an optimum model SED curve, for the redshift and 
age of the stellar population simultaneously, was car- 
ried out, and the derived values were compared. 

Thus, we estimated both the galaxy age and the 
redshift within the frames of the specified models (see 
Verkhodanov et al., 1998a, 1999). It is clear from gen- 
eral considerations that the reliability of the result at 
large redshifts shows a significant dependence on the 
availability of infrared data (up to the R range), be- 
cause when fitting, we "cover" the region of rapid 
spectrum change (a jump) before the optical range 



of SED and thereby can reliably (with a well de- 
fined maximum of the plausibility curve) determine 
the location of our data. This can be proved as fol- 
lows. When the available points are removed (to check 
the reliability of the procedure) with retaining only 
3 points (one of which is in the K range), we obtain 
in fitting the same result on the curve of deviations 
as for 4 or 5 points. If the infrared range is not used, 
the result will then be more uncertain. However, as 
we show (Verkhodanov et al., 1999) the version of 
"tight" positioning of the four filters, as in our case 
of BVRI photometry, gave a good result in a sample 
of 6 galaxies, which is a good fit to the result obtained 
with using all the filters, including the infrared range. 

A similar procedure was employed to estimate the 
rcdshift and the age for galaxies of the RC catalogue. 
The common distinction from the papers by Verkho- 
danov et al., 1998b, 1999) is that our observations 
were made only in 4 filters, BVRI, so, the infrared 
data are lacking. 

Table 4 presents redshift and age estimates for the 
stellar systems of radio galaxies from the RC cata- 
logue. The rcdshift z c , the age of the radio galaxy 
t, the discrepancy e and the age of the Universe T 
at z — z c (H = 65, il m — 0.3, fl\ = 0.7 for stan- 
dard Friedmann-Lemaitre- Robertson- Walker cosmol- 
ogy (Thomas & Kantowski, 2000)) are given in the 
table columns for the models GISSEL and PEGASE, 
respectively. 



Table 4: Colour redshifts and ages of stellar systems for the RC Catalogue radio galaxies 



Object RC J 


GISSEL'98 


PEGASE 


Zc 


t[Myr] 


e 


T[Myr] 


Zc 


t[Myr] 


£ 


T[Myr] 


1 


2 


3 


4 


5 


6 


7 


8 


9 


0015+0503a 


0.73 


900 


0.0412 


7546 


0.48 


3000 


0.0299 


9237 




0.07 


16000 


0.0342 


13512 


4.91 


16000 


0.0199 


1273 


0015+0501 


0.81 


1000 


0.0334 


7103 


0.87 


3250 


0.0358 


6797 


0034+0513 


0.98 


16000 


0.0729 


6286 


1.29 


10000 


0.0374 


5130 












1.11 


15000 


0.0196 


5755 


0039+0454 


0.95 


16000 


0.0346 


6419 


0.99 


7000 


0.0057 


6242 


0105+0501 


0.28 


500 


0.0852 


11032 


4.21 


9250 


0.0839 


1537 


0135+0450 


0.35 


2500 


0.0106 


10350 


0.10 


7250 


0.0018 


13112 




0.12 


16000 


0.0072 


12855 










0152+0453 


0.75 


400 


0.0051 


7431 


0.81 


640 


0.0016 


7103 




1.15 


900 


0.0113 


5606 










0159+0448 


0.13 


2000 


0.0093 


12729 


0.41 


1900 


0.0093 


9813 












0.09 


5000 


0.0090 


13244 












4.75 


15000 


0.0060 


1326 


0209+0501a 


0.39 


1800 


0.0067 


9988 


0.39 


3500 


0.0063 


9988 












0.09 


6750 


0.0088 


13244 












0.07 


9250 


0.0097 


13512 
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Table 4: Colour redshifts and ages of stellar systems for the RC Catalogue radio galaxies (continued) 



1 


2 


3 


4 


5 


6 


7 


8 


9 














0.07 


16000 


o 


0058 


13512 


0209+0601b 


0.69 


13000 


o 


0400 


7783 


0.70 


6000 


o 


0483 


7723 


0318+0456 


0.67 


16000 


n 

u 


0491 


7906 


0.71 


6500 


o 

U 


0071 


7663 


0444+0501 


1.17 


200 


n 

u 


0365 


5534 


2.42 


13000 


o 

u 


041 9 

U4:l iJ 


2870 


0457-1-0452 


0.40 


2000 


n 

u 


0092 


9900 


0.41 


3500 


o 


0068 
uuuo 


9813 


0836+051 1 


0.77 


250 


n 

u 


0058 


7319 


0.81 


286 


o 

U 


0094 


7103 




1.14 


600 





0028 


5643 


1.21 


2100 


o 


0037 


5394 




1.93 


2500 





0229 


3599 


1.37 


2400 


o 


0021 


4886 


0837+0446 

V / V < 7 f | V I 1 1 V 1 


0.98 


200 


n 

u 


0490 


6286 


0.99 


202 


o 

VJ 


0494 


6242 




1.91 


6000 





0697 


3636 


1.92 


9250 


o 


0635 


3617 


0845+044 

V I V • I \7 | V 1 -L -L 


5.29 


16000 


n 

u 


1 1 39 


1160 


0.66 


8750 


o 

VJ 


0560 

U<JUU 


7969 


0908+0451 

Ul/UU 1 U ill 1 


0.49 


1400 


n 

u 


0322 


9159 


0.48 


3250 


o 

VJ 


01 93 


9237 














4.99 


16000 


o 


0125 


1248 


0909+0445 

V 1 - 1 V 7 - 1 | V I 1 1 1 


0.67 


1800 


n 

u 


01 45 

\J J_4:U 


7906 


0.64 


4000 


n 


01 fiS 

U1UO 


8096 


0934+0505 

VI - 1 '.7 1 | V7*^v7*-/ 


1.72 


1200 


n 

u 


0785 


4010 


1.79 


4250 


n 

VI 


0728 


3865 














1.79 


16000 


o 

VJ 


0727 

U 1 J 


3865 


1011+0502 

1 V 1 1 | V7 V7 v7 


0.46 


300 


n 

u 


0^01 

\JO\J ± 


9397 


0.50 


453 


n 

VJ 


0624 


9082 


1031+0443 


0.86 


2500 


o 

u 


01 88 


6846 


0.87 


4250 


o 

VJ 


UZj J o 


6797 


1043+0443 


1.19 


3000 


n 

u 


01 38 


5463 


0.72 


2200 


VJ 


UZj<JU 


7604 














1.23 


4500 


o 

VJ 


01 71 

Ull 1 


5326 














5.98 


8500 


n 

VJ 


01 81 


993 


1124+0456 

i i z_ i | v i l v i V/ 


0.35 


10000 


n 

u 




10350 


0.35 


5750 


n 


U -L -L <J 


10350 


1142+0455 

1 1 1 Z_ | V 1 1 '.. 1 '.. 1 


5.20 


16000 


n 

u 


1 082 


1185 


4.99 


571 


n 

VJ 


uoou 


1248 














0.35 


14000 


VJ 


0421 


10350 


1152+0449 

1 17 — | V l 1 1 > 7 


1.17 


16000 


n 
u 


Ofil Q 


5534 


1.33 


8000 


n 

VJ 


uuuy 


5006 














1.30 


15250 


n 

VJ 


D07Q 

UU 1 t7 


5099 


1155+0444 

7 1 U U | v7 J. J. J. 


4.18 


500 


n 

u 


1161 


1550 


0.33 


5750 


o 

VJ 


01 79 


10538 




0.35 


9000 


n 

u 


0091 


10350 












1213+0500 

ItilU I V7 V7 V7 V7 


0.70 


700 


n 

u 


0201 


7723 


0.68 


2100 


o 

VJ 


U 1UZ; 


7844 




1.08 


2000 


n 

u 


0334 


5871 


5.15 


15500 


VJ 


01 1 8 


1200 




0.07 


9000 


n 

u 


0522 


13512 












1235+0435b 

vj | V7 j: V7 V7 U 


0.64 


11000 


n 

u 


0042 


8096 


0.68 


5750 


o 

VJ 


0083 
uuoo 


7844 


1322+0449 

LfJ Z_' Z_ | V I 1 1 - 7 


0.65 


16000 


n 

u 


1129 


8032 


0.78 


7250 


o 

VJ 


0702 

U 1 UZj 


7264 














0.99 


14000 


o 

VJ 


0037 

UUO 1 


6242 


1333+0452 


1.04 


1400 


n 

u 


01 58 


6032 


1.08 


4000 


o 

VJ 


0062 

UUUZj 


5871 


1339+0445 

J < 7 '.7 - 7 | V I 1 1 1 


0.67 


4500 


n 

u 


021 7 


7906 


0.71 


5000 


o 

VJ 


01 36 

U -LOU 


7663 


1357+0453 

1 V7 .7 1 J V7 j: V7 V7 


0.73 


1400 


n 

u 


0421 


7546 


0.41 


4000 


VJ 


0409 


9813 


1429+0501 


5.38 


450 





0827 


1135 


0.61 


11500 





0483 


8294 


1436+0501 


1.35 


800 





0157 


4945 


1.15 


1680 





0136 


5606 




1.93 


3000 





0318 


3599 


1.18 


1900 





0127 


5499 


1446+0507 


0.35 


5000 





0140 


10350 


0.24 


5750 





0191 


11451 














4.99 


6750 





0189 


1248 














0.11 


14500 





0144 


12983 


1503+0456 


0.78 


500 





0065 


7264 


0.83 


1015 





0023 


6998 


1510+0438 


0.66 


8000 





0110 


7969 


0.75 


5750 





0074 


7431 


1551+0458 


1.05 


1200 





0341 


5991 


1.11 


4000 





0294 


5755 




1.32 


6000 





0099 


5036 


1.32 


5750 





0129 


5036 


1626+0448 


0.04 


200 





1268 


13930 


0.03 


202 





1282 


14073 




2.31 


6000 





1642 


3011 


2.30 


3000 





1814 


3025 
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Table 4: Colour redshifts and ages of stellar systems for the RC Catalogue radio galaxies (continued) 



1 


2 


3 


4 


5 


6 


7 


8 


9 


1638+0450 


0.84 


500 


n 


031 5 

UOl o 


6947 


0.89 


1015 


n 

VJ 




6699 




1.67 


2000 


o 


0284 


4119 


1.78 


4750 


o 


0134 


3885 


1646+0501 


0.65 


>5000 


n 


1 1 71 

1111 


8032 


0.64 


6750 


n 


01 06 

VJ 1UU 


8096 


1703+0502 


0.81 


600 


n 

VJ 


0692 


7103 


3.56 


640 


n 


0567 


1874 




0.07 


2500 


n 

\J 


0431 


13512 


0.35 


1680 


n 

\J 


0508 

VJOVJO 


10350 




0.00 


6000 





0452 


14515 


0.00 


5500 





0363 


14515 


1706+0502 


0.96 


1800 





0339 


6374 


1.06 


5750 





0293 


5951 




1.17 


16000 





0317 


5534 












1722+0442 


1.02 


4000 





0033 


6115 


1.01 


4500 





0033 


6157 


2029+0456 


0.81 


800 





0094 


7103 


0.88 


2600 





0107 


6747 


2219+0458 


0.98 


1200 





1271 


6286 


1.36 


7750 





0555 


4916 




1.22 


16000 





0987 


5360 












2224+0513 


0.73 


1200 





0101 


7546 


0.77 


3500 





0072 


7319 


2247+0507 


0.78 


700 





0094 


7264 


0.93 


2500 





0068 


6510 




1.38 


3000 





0253 


4857 


1.45 


4750 





0134 


4660 


2348+0507 


1.43 


450 





0016 


4715 


1.24 


571 





0030 


5293 














2.07 


3000 





0291 


3362 



Notes to Table 4 

0034+0513 — the age of the stellar system in the 
two alternatives from the model GISSEL'98 is greater 
than that of the Universe at the given z; we choose the 
age where the difference is smaller despite a somewhat 
greater error. 

0039+0454 the age of the stellar system 
from the PEGASE model is inadmissivcly large, it 
is about 3 times as large as that of the Universe 
at a given z. The same is observed in a number of 
other sources (0209+0501, 0845+0454, 1142+0455, 
1152+0449, 1235+0435b, 1322+0449, 1510+0438, 
1646+0501, 2219+0454). The cause of this has not 
been understood yet. 

0105+0501 in contrast to the majority of ob- 
jects of our sample the redshift and the age from the 
two models are not consistent. Besides, the computa- 
tion gives large errors in this case. The optical images 
(Soboleva et al. , 2000) show that the bright line Lya 
falls in the V filter. For this reason, the stellar mag- 
nitude is strongly corrupted. 

0135+0450 — the redshift value lies between 
0.34 and 0.10, but closer to 0.34. Note that from 
the GISSEL'98 model there is another minimum at 
z = 0.17 with the age 5750 Myr and a deviation of 
0.0105. The values at the absolute minimum are dis- 
carded in the two models because of the discrepancy 
between the age of the system and that of the Uni- 
verse. 

0159+0448 — z = 0.09 (GISSEL'98) and z = 
0.13 (PEGASE) lie outside the region of permissible 
values in the plane z — R (see Fig. 1; Pursimo et al., 
1999). The age 15000 Myr in the GISSEL model does 
not conform to reality. We consider the alternative 
with z = 0.41 to be more reliable. 



0209+0501a z = 0.39 for both models. The 
age derived does not conform to real fact for a min- 
imum deviation 0.0058 in the GISSEL model. Two 
other estimates yield a larger value of deviation. 

0318+0456 — close redshifts for the models GIS- 
SEL'98 (z = 0.71) and PEGASE (z = 0.67). The age 
turns out to be ultimately small for both models. 

0444+0501 — although to the value z = 1.17 
for the PEGASE model corresponds a minimum de- 
viation, the closest to the spectroscopic redshift (z = 
2.73) is obtained in the GISSEL model (z = 2.42). 
The estimated age proves to be too great to be real, 
although z is close to real. 

0836+0511 — the values z = 1.21 (GISSEL'98) 
and z = 1.14 (PEGASE) are close (to an accuracy of 
6 %) and practically have a minimal deviation. One 
can see here that the age of stellar systems turns out 
to be older in the GISSEL models. 

0845+0444 — z = 5.29 (PEGASE) is not admis- 
sible both by age and by position in the plane z — R. 

0908+0451 — the colour redshift (z ph = 0.48) 
is confirmed by spectral measurements at the 6 m 
telescope of SAO RAS (z sp = 0.5). 

1142+9455 — z = 4.99 (GISSEL'98) and z = 
5.20 (PEGASE) do not fall within the range of ad- 
missible values in the plane z— R. We choose z = 0.35, 
though in this case the age of the stellar system, ac- 
cording to the criterion of minimum of squares of de- 
viations, proves to be somewhat larger than the age 
of the Universe. 

1213+0500 — z = 5.15 (GISSEL) is outside of 
the region of acceptable values in the plane z — R, and 
the age of the stellar system in this case is 10 times 
that of the Universe. There are close redshift values 
in the region z — 0.7 for both models, a minimum x 2 



RESULTS OF INVESTIGATION OF RADIO GALAXIES OF THE SURVEY "COLD 



13 



being consistent with this solution in the PEGASE 
model. 

1322+0449 — the variant with z = 0.99 in the 
GISSEL model is unsuitable because the age of the 
stellar system is twice that of the Universe. A similar 
situation also occurs in the PEGASE model where 
the solution goes beyond the limits of acceptable age 
values. There is a stable minimum for z = 0.78 in the 
PEGASE model. In all the cases the redshift proves 
to be less than 1. 

1357+0453 — assume z = 0.73 from PEGASE 
model, z = 0.41, however, for GISSEL model falls out 
of the plane z — R. Note that there is one more local 
minimum for the GISSEL model, which corresponds 
to z = 3.86 in the plane of distribution of deviations, 
but the value of this minimum is 3 times as large as 
for z = 0.41, and the corresponding point also falls 
out of the plane z — R. 

1429+0501 — we fail to find acceptable values of 
the parameters for either model. This may be due to 
difficulties of separating a radio galaxy from a nearby 
star in the optical range (Parijskij et al., 1996a). 

1436+0501 — accept the alternatives z = 1.18 
(GISSEL) and z = 1.35 (PEGASE) with errors close 
to minimum. The difference between the redshifts in 
these models is then about 12 %. 

1446+0507 — at z = 0.11 the age of the galaxy 
in the GISSEL model turns out to be larger than the 
age of the Universe at a minimum of deviations, at 
z = 4.99 with a minimum error it does not fall within 
the range of admissible values in the plane z — R for 
radio galaxies, and the age of the stellar system is 5 
times the age of the Universe. We retain z — 0.35 
for PEGASE and z = 0.24 for GISSEL as possible 
alternatives. 

1626+0448 — spectroscopic redshift for this ob- 
ject is z = 2.66 (Afanasiev et al., 2002). The cur- 
rent calculations have been carried out with no pro- 
vision for the line Ly a. The data obtained show that 
z = 0.03 and z = 0.04 do not fall in the interval of 
acceptable values for radio galaxies in the plane z — R. 
In all the cases the large errors are due to the fact that 
the bright line Ly a falls on the B band (Parijskij et 
al., 1996a). The z values equal to 2.31 for PEGASE 
and 2.30 for GISSEL are quite close to the true val- 
ues. When the line Ly a is taken into account in the 
computations, the redshift value then turns out to be 
equal to the spectroscopic. 

1638+0450 — choose z = 1.78 (GISSEL'98) and 
z = 1.67 (PEGASE), though z = 0.89 (GISSEL'98) 
and z = 0.84 (PEGASE) cannot be rejected either. 

1703+0502 — the stellar magnitudes in all the 
filters can be affected by a nearby bright star (Pari- 
jskij et al., 1996a). 

1722+0442 the spectroscopic z = 0.7 
(Afanasiev et al., 2002) is by 30% lower than the 
value derived by photometry for both models. 
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15 
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1.1 



1.2 
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Figure 2: Histogram characterizing the distribution 
of the ratio of colour redshifts from the GISSEL '98 
(zgi SS ) and PEGASE (z peg ) models. The central part 
of the histogram, within which 80 % of objects fall, is 
presented. 



2247+0507 — there are two minima in the two 
models. The z values less than 1 differ for the models 
by less than 20%. The difference is less than 10% for 
z of about 1.4, however, the age of the stellar system 
is the same as the age of the Universe for the model 
GISSEL. By the criterion x 2 we choose the former 
alternative in both models. 

2348+0507 — the stellar magnitudes in all the 
filters can be affected by a nearby bright star (Pari- 
jskij et al., 1996a). 

6. Discussion and results 

1. The problem of colour redshifts of distant 
galaxies 

As it was shown earlier (Verkhodanov et al., 
1999), the model PEGASE yielded satisfactory re- 
sults over 40 objects for which spectral measurements 
of redshifts are available. Analysis of Table 4 shows 
that the new model CISSEL often gives results only 
slightly differing from the PEGASE model (see also 
Verkhodanov et al., 2001a, b). 

Fig. 2 presents a histogram of the ratio (1 + 
zgissel)/{1 + zpegase)- The distribution of errors 
is not normal, there is a nucleus with very small model 
errors, in which most of the objects (80%) fall, but 
20% of the objects have very large errors. As a rule, 
the large errors are caused by getting of the strong 
line {Ly a) in one of the filters, illumination from a 
nearby bright object and complicated cases of SED 
deviation from the model. A number of such exam- 
ples are described in notes to Table 4. 
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It should be noted that we use the standard set 
of filters of the 6 m telescope, BVRI, which is sat- 
isfactory for not too distant objects (the colours of 
which have been measured yet). Measurements have 
to be extended to the H and K ranges. In a consid- 
erable number of cases, this will enable elimination 
of uncertainties in the estimates. The situation may 
change for ultimately distant objects with R> 24 m , 
and we have to wait until the observations are com- 
pleted, prior to giving recommendations. It is clear 
from general considerations that there is a danger 
from "right" and "left", that is, the secondary star 
formation sites may distort the blue region of the 
spectrum, but the dust at very great redshifts may 
deform the IR region. It can easily be shown that 
large discrepancies may point to possibilities of lu- 
minosity distortion by strong lines, and we hope to 
take this into account in the future. For z > 2, it is 
necessary to take into account, at least, the line Ly a. 
Considerable errors may arise because of the periodic- 
ity of series in Bor's model, but here simple energetic 
considerations may be helpful (discrepancies in the 
colour and photometric estimates). This can be illus- 
trated by the RC object J1703+0502 formally having 
a zero redshift. When adopting this estimate, the op- 
tical luminosity of this object will then prove to be 
so low that it cannot produce any perceptible radio 
radiation (note that P ra< n ^L^ t ; Iskudarian & Par- 
ijskij, 1964; Franceschini et al., 1998). When making 
use of the secondary criteria, we sometimes have to 
reject the variant with the smallest discrepancies and 
take the next one. 

2. The age of stellar systems of host galaxies 
As it is known, the age is much more difficult to 
assess than the redshift. The older the stellar popu- 
lation the larger the error may be. We have made a 
histogram (Fig. 3) of differences of ages determined 
from the models GISSEL (t glss ) and PEGASE (t peg ) 
for one and the same object. They are normalized 
to the age of the Universe with A = 0.7 (T) for the 
moment that corresponds to the measured redshift 
(tgi S s — t peg )/2T. The distribution of differences of 
ages derived from different models are far from nor- 
mal, but there is a nucleus where 70 % of objects fall. 
The "model" dispersion for them is close to 20%. 
Large departures were detected in 9 out of 50 ob- 
jects with the use of PEGASE. The measured age of 
host stellar systems for them is far larger than that 
of the Universe even at redshifts close to GISSEL. 
For this reason, we consider the GISSEL model to 
be more trustworthy. A systematic difference between 
the models has been discovered, which increases with 
rising redshift (see Fig. 4). When the difference is 
taken into account, the dispersion of ages drops to 
< 10 % for objects with close redshifts in both mod- 
els. 

It is essential that not a single GISSEL model 



N 
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5- 
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Figure 3: Histogram of half- difference of colour ages 
of host stellar systems obtained from the models GIS- 
SEL'98 {tgiss) and PEGASE (t peg ) vrs the age of the 
Universe with A = 0.7 (T) for the moment corre- 
sponding to the measured redshift. The central part of 
the histogram is shown where 80 % of objects fall. 

galaxy is older within the errors than the Universe in 
the model with the A-term equal to 0.7-0.8. There are 
more than ten such galaxies in the old model SCDM, 
and they can be used for independent assessment of 
the behaviour of H(z), A-term and quintessence (sug- 
gestion of Starobinsky (1998), see Saini et al. (2000); 
Parijskij et al. (1998)). The average age of galaxies is 
a few billion years, however, there is a subgroup of 
galaxies having the age that coincides within the er- 
rors with the age of the Universe at a corresponding 
redshift, and a group of formally very young (some 
hundred million years) stellar systems. The formers 
are likely to belong to the first generation of galaxies 
in the Universe which formed at z > 1 (z ~ 5 — 10). 
The situation of the galaxies with "young" objects is 
more complicated, their colour can be distorted by 
repeated star-formation bursts at merging of galax- 
ies or under the action of close passings. The variant 
of "young" galaxy cannot be accepted, since we deal 
with powerful radio galaxies, for functioning of which 
a supermassive black hole is necessary (DMO — Dark 
Massive Object, Salucci et al., 1999) with the mass 
of about 10 9 M Q , that is not possible within the stan- 
dard models of black hole formation (Franceschini et 
al., 1998). Only the primary black holes with masses 
10 4 — 1O 6 M , round of which galaxies form later on, 
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Figure 4: Systematic difference between ages of stellar 
systems from GISSEL'98 (t glss ) and PEGASE (t peg ) 
models normalized to the age of the Universe with 
A = 0.7 (T) for the moment of the corresponding 
redshift as a function of redshift. The departures from 
the relationship derived are 6 %. 



may be the alternatives of "merging" . 

As a result, a careful conclusion can be drawn that 
at least the statistical estimates of the redshift and 
age for population of powerful radio galaxies give sat- 
isfactory results. The GISSEL model can be recom- 
mended not only for radio-quiet galaxies but also for 
powerful radio galaxies. It is advantageous to use all 
available data on the objects to obtain more reliable 
estimates, this will decrease the errors. 

In conclusion we enumerate the results obtained. 

1. The largest currently available body of data 
on BVRI magnitudes of host galaxies responsible for 
the origin of powerful radio galaxies is presented. 

2. The colour redshifts for powerful radio galax- 
ies show a satisfactory agreement with the spectral 
ones (the error is 10-20% with a small fraction of 
large errors). The recent 2001 spectral observations 
of the RC objects RC J0908+0451, RC J1154+0431, 
RC J1626+0448 RC J1722+0442 gave errors of red- 
shift measurement by the techniques described above 
within 10-15% (Afanasiev et al., 2002). 

3. The limited number of closely spaced filters as 
in our BVRI case may also yield satisfactory results 
even for large redshifts. 

4. The redshift distribution for the studied ob- 
jects of our sample (the subgroup of objects brighter 
than mji = 23. 5 m ) shows a maximum near z~l, i.e. 
in the range of maximum radio activity of the Uni- 
verse. The group of objects with large colour redshifts 
(z > 2.5 — 3) requires a separate analysis. In any case, 
we do not consider yet the gap in the population of 
the region (1.5 < z < 2.5) to be real. 

5. The colour data are generally not at variance 
with the stellar magnitudes in the filter R (Parijskij 



et al., 1996a) when R < 22. 5 m . Fainter objects show 
a higher dispersion of photometric redshifts at one 
and the same value of R. One can notice here two 
branches in the plane (z — R) (Pursimo et al., 1999). 
The search for differences in the morphology, radio 
luminosity, spectral indices has given no final result. 
Note that the objects with the steepest spectra and 
a large radio-to-optical luminosity ratio occur in the 
branch with large redshifts. However, further inves- 
tigations are needed. Objects with low relative radio 
luminosity, as it was to be expected, prove to be either 
quasars or nearby galaxies. 

6. The age of galaxies is determined less reliably, 
and results turn out to be of low significance for larger 
z. However, in practice one can always indicate the 
lower limit of the age of galaxies and, therefore, the 
minimum redshift of their formation. This age is al- 
ways larger than the standard estimate of the object 
lifetime, and in a number of cases it exceeds the age 
of the Universe at an object redshift in a simple CDM 
model. There are no such galaxies in a model with a 
A-term of 0.6-0.8 (e.g. de Bernardis et al., 2000). 

7. A programme has been developed and tested 
of automatic determination of colour redshifts and 
ages of galaxies, which is applicable to any redshifts 
with provision for transformation of the shape of the 
filters when changing over from the rest to the mov- 
ing system of reading, available through the server 
http://sed.sao.ru . 



We contemplate further investigation of colour and 
photometric techniques as applied to the population 
of distant radio galaxies. New more realistic models 
of colour evolution and refined methods of age evalu- 
ation in stellar systems will make it possible to obtain 
more reliable results for a great number of objects. 
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Fig. 5. Oprical identifications of RC sources: VLA radio isophotes overlaid on the 6 m telescope R-band 
images, and R-band isolines overalaid on the VLA images. 
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Fig. 6. Fitted SED models and corresponding probability functions for RC objects. 
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